Coupled multimode optomechanics in the microwave regime 
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The motion of micro- and nanomechanical resonators can be coupled to electromagnetic fields. 
This allows to explore the mutual interaction and introduces new means to manipulate and control 
both light and mechanical motion. Such optomechanical systems have recently been implemented 
in nanoelectromechanical systems involving a nanomechanical beam coupled to a superconduct- 
ing microwave resonator. Here, we propose optomechanical systems that involve multiple, coupled 
microwave resonators. In contrast to similar systems in the optical realm, the coupling frequency 
governing photon exchange between microwave modes is naturally comparable to typical mechanical 
frequencies. For instance this enables new ways to manipulate the microwave field, such as mechan- 
ically driving coherent photon dynamics between different modes. In particular we investigate two 
setups where the electromagnetic field is coupled either linearly or quadratically to the displacement 
of a nanomechanical beam. The latter scheme allows to perform QND Fock state detection. For 
experimentally realistic parameters we predict the possibility to measure an individual quantum 
jump from the mechanical ground state to the first excited state. 

PACS numbers: 85.85. +j, 84.40.Dc, 42.50.Dv 



Introduction. - Significant interest in the interaction 
and dynamics of systems comprising micro- and nanome- 
chanical resonators coupled to electromagnetic fields, as 
well as the prospect to eventually measure and control 
the quantum regime of mechanical motion, has stimu- 
lated the rapidly evolving field of optomechanics (see [1] 
for a recent review). In the standard setup, the light 
field, stored inside an optical cavity, exerts a radiation 
pressure force on a movable end-mirror whose motion 
changes the cavity frequency and thus acts back on the 
photon dynamics. This way, the photon number inside 
the optical mode is linearly coupled to the displacement 
of a mechanical object. Beyond the standard approach, 
new developments have introduced optical setups with 
multiple coupled light and vibrational modes pointing 
the way towards integrated optomechanical circuits [2]- 
[6]. These systems allow to study elaborate interactions 
between mechanical motion and light such as mechani- 
cally driven coherent photon dynamics that introduces 
the whole realm of driven two- and multi-level systems 
to the field of optomechanics [7]. Coupled multimode 
setups furthermore allow to increase measurement sen- 
sitivity j8j and enable fundamentally different coupling 
schemes. Accordingly, recent experiments achieved cou- 
pling the photon number to the square and quadruple of 
mechanical displacement [2", 9 . Such different coupling 
schemes are needed, for instance, to afford quantum non- 
demolition (QND) Fock state detection of a mechanical 
resonator [2l,[T0Hl2]. 

Besides optics, recent progress has made it possible to 
realize optomechanical systems in the microwave regime 
[I3l m] . In this case the optical cavity is replaced by a 
superconducting microwave resonator whose central con- 



ductor capacitively couples to the motion of a nanome- 
chanical beam. This optomechanical approach consti- 
tutes a new path to perform on-chip experiments mea- 
suring and manipulating nanomechanical motion that 
adds to electrical concepts using single electron transis- 
tors [TSHiT], superconducting quantum interference de- 
vices [ISl [19] , driven RF circuits [20] or a Cooper-pair- 
box [21] [22] . One advantage of on-chip optomechanics is 
to use standard bulk refrigerator techniques in addition 
to laser cooling schemes [23j [24] . This recently enabled 
cooling a single vibrational mode close to the quantum 
mechanical ground state [25]. Furthermore, nonlinear 
circuit elements can be integrated. This afforded ultra- 
sensitive displacement measurements with measurement 
imprecision below that at the standard quantum limit 

Here we go beyond single mode systems, that have 
been considered for optomechanics in the microwave 
regime so far, and propose setups with coupled mi- 
crowave resonators. For these systems, the coupling fre- 
quency between microwave modes turns out to be com- 
parable to the mechanical frequency. This has several 
implications both for the classical and the prospective 
quantum regime. As an example, focusing on the cur- 
rently in experiments accessible regime of classical mo- 
tion, we demonstrate how this allows to manipulate the 
microwave field in terms of mechanical driving. 

Two-resonators setup with linear mechanical 
coupling. - We consider the coplanar device geome- 
tries depicted in Fig. [l^ with two identical superconduct- 
ing microwave resonators a/,, (ir. The central conduc- 
tors of aL and aR are assumed to adjoin for a length dg 
that is much smaller than the total wave guides' length 




FIG. 1. Schematic device geometry for two superconducting 
microwave resonators aL, clr with a nanomechanical beam 
coupled to aL ■ (a) The two resonators (each of length d) are 
coupled to external feed and transmission lines (green), (b) 
The central conductors of ql and an (red) capacitively couple 
due to a small region of length dg where the resonators adjoin. 

(c) At the other end of aL a small mechanical beam, connected 
to ground (blue), is placed. Its displacement x affects the line 
capacitance of resonator aL changing its resonance frequency. 

(d) System's resonance frequency as function of displacement: 
the beam's displacement x linearly changes the bare mode fre- 
quency of aL while the one of an is unaffected (dashed). Due 
to the coupling g between modes there is an avoided cross- 
ing 2g in the eigenfrequencies (blue), (e) Analogous optical 
setup: a static, dielectic membrane, placed inside a cavity 
with a movable mirror, couples two separate optical modes 
aL, aR. 



d [Fig. fit)]. A nanomechanical beam, connected to the 
ground plane, is placed at the other end of aL [Fig. Ilfc]. 
Its motion in terms of displacement x changes the line 
capacitance c (capacitance per unit length) between the 
central conductor and the ground plane in a small region. 
In the following, we will derive and discuss the Hamil- 
tonian for the system depicted in Fig. [l] starting from 
a single microwave resonator whose line capacitance is 
changed due to the motion of a mechanical beam [Fig. lip] 

Single MW resonator coupled to a nanomechanical 
beam. - We concentrate on one of the microwave res- 
onators in Fig. [l] The Lagrangian of an electric circuit 
can be conveniently expressed in terms of a flux variable 
0(x, t) = J_ dr V{x^ r), where V{x^ t) = dt(l){x^ t) is the 
voltage on the transmission line at position x and time 
t, see for instance [27 . For a finite length microwave 
resonator with line inductance / and line capacitance c 
(both per unit length) , the corresponding Euler-Lagrange 
equation yields a wave equation with speed v = y/l/lc. 
Using appropriate boundary conditions, ^(x, t) can be ex- 



panded into normal modes 0^, see [28]. The Lagrangian 

reads C = ^^/^[f^^ + ^ (^) ^k\ ^^^ ^^^ microwave 
resonator possesses resonance frequencies uj^ — kiiv/d. 
The Hamiltonian is obtained by Legendre transformation 
using the canonically conjugated momentum ti^ = c^k- 
Quantizing the system by introducing creation and anni- 
hilation operators a|., ak for the individual modes /c, with 
[a/c, aj./] = Sk^k'j the Hamiltonian is a sum of harmonic 
oscillators, ^^ fyjJkCi\ak^ where we neglected the vacuum 
energy. 

For the resonator aL in Fig. [T] the motion of the me- 
chanical beam will change the line capacitance c at the 
end of the wave guide. If we denote the total capacitance 
between the central conductor and the beam by Ch and 
define the optomechanical frequency pull per displace- 
ment, e/c = —duok/dx^ we find 

^'^^Z-a;^/(2^fe), (1) 



e/c 



dx 



where we used the line impedance Z = \/l/c. Up to 
linear order in x, the Hamiltonian for a single microwave 
resonator coupled to a nanomechanical beam reads H = 
^f,h{uk,o -ekx)alak. 

Coupled MW resonators. - We now turn to the cou- 
pled resonator setup in Fig. fl] As the length dg of the 
region where aL and aR adjoin is much smaller than 
the resonators' length d, the capacitive coupling between 
both central conductors can be considered in terms of 
a constant, total capacitance G. We briefly switch to 
a discretized description. If we neglect the motion of 
the nanomechanical beam the circuit diagram looks as 
depicted in Fig. [2J where the flux variable ^(x, t) is dis- 
cretized as (pnit) at position x = na, with a being the con- 
stant spacing between nodes. For each individual node, 
an equation of motion can be written down [27 . From 
this, taking into account appropriate boundary condi- 
tions for a ^ 0, the continuous version of the correspond- 
ing Lagrangian is found to be 
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where ^l,/c, 4^R,k refer to the individual normal modes of 
the left and right resonator, respectively. The first two 
terms describe two separate wave guides, while the last 
term characterizes the coupling between both. 

To transform to the Hamiltonian, we consider 
the canonically conjugated momentum 7rL[R]^k = 

dL/d(j)L[R],k = C(j)L[R],k T '^^{T.k^R.k - T.k^L,k)' In 
the following we will restrict to a single mode in each 
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FIG. 2. Discretized circuit diagram for the system depicted 
in Fig. [l] without nanomechanical beam. The two central 
conductors are coupled via a total capacitance G. I and c 
are the line inductance and capacitance, respectively, a is the 
spacing between discrete points. 



resonator {k = 1), and drop the label referring to the 
mode index. For G/d <C c (see discussion below), 
we can simplify the expression for tt^ji^^j. and consider 

^L[R],k = C(j)L[R],k' The first two terms of J2| transform 
into two harmonic oscillators of frequency ul and ujr. 
For the coupling we have to consider {^r — ^l)'^ with 
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Using rotating wave approximation, we find the Hamil- 
tonian, 



H = ncjL 1 
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where we neglected the vacuum energy. Note that none 
of the resonators in Fig. [l] is short-circuited such that 
there are voltage antinodes at both ends of each resonator 
allowing to have maximal coupling between a^ and aR^ 
as well as to the feed and transmission line. 

The coupling between the resonators has two effects. 
First, both frequencies cjl, ^r^ originally defined for un- 
coupled modes, are lowered by a constant value. In the 
following, this shift of frequency is neglected by simply 
redefining the resonators' frequencies. More important 
is the coupling between modes in terms of the coupling 
frequency 



g = ^ulujr — . 



(3) 



Finally we take into account the motion of the nanome- 
chanical beam changing the left mode's bare eigenfre- 
quency uoq in the way discussed above, uol{x) = uoq — ex 
(see Eq. (IT]) with /c = 1). The final Hamiltonian for the 
system depicted in Fig. [T] then reads 



% = h (cjo — ^x) o/^aL + huJoa^j^aR 
+ hg [a^^aR 



CLrCIL 



(4) 



In principle, according to Q with uol{x)^ the coupling 
frequency g depends on displacement x. However, for 
typical parameters the dependence is negligible and g can 
be considered to be constant. The resonance frequency 
of Q is depicted in Fig. [l]i. 



Coupling frequency comparable to the mechanical fre- 
quency (g ^ rt). - Given Eq. (|3|, the coupling fre- 
quency between the two resonator modes reads g/ujQ = 
{cg/c) • (dg/d)^ where we defined the coupling line ca- 
pacitance Cg = G/dg along the length of the coupling 
region dg. In general, Cg will be much smaller than 
the line capacitance between each central conductor and 
the ground plane c: first of all, the distance between 
the two central conductors is significantly larger than 
the distance between a single conductor and the adja- 
cent ground plane. Second, the capacitance between 
the central strip lines is shielded by the grounded re- 
gion in between. Here we crudely assume Cg/c = 10~^. 
For d in the cm range and dg c^ 0.1mm {dg is chosen 
such that a several 10 /im long nanomechanical beam 
can be fabricated in between the region where the res- 
onators align), we have dg/d = 10~^ and the coupling 
between modes is g/ujo = 10~^ where uq will be in the 
GHz range. Common eigenfrequencies of nanomechan- 
ical beams are Q = 100 kHz - 10 MHz. Hence, due to 
their much smaller photon frequency, coupled multimode 
optomechanical systems in the microwave regime natu- 
rally possess coupling frequencies in the range of typical 
mechanical frequencies {g c^ Q). The relevance of this 
regime for instance to realize all kinds of driven two- and 
multi-level photon dynamics in optomechanical systems 
has been pointed out in [7 . 

Transmission spectrum. - As an example to emphasize 
the characteristics of coupled optomechanical systems in 
the microwave regime and to demonstrate implications of 
g c^ Q even in the presently accessible regime of classic 
mechanical motion, we will discuss how the microwave 
field in the setup of Fig. [l] can be manipulated in terms 
of mechanical driving (see [29] for a universal mechanical 
actuation scheme). Experimentally, the impact can be 
most easily observed in terms of the transmission spec- 
trum. We assume the left resonator a^ to be driven at 
frequency ujl via the feed line, while the transmission 
down the transmission line is recorded. We consider the 
coupling of the left (right) resonator to the feed (trans- 
mission) line in terms of the the resonators' decay rate K:. 
Given the Hamiltonian Q, using input /output theory, 
the equation of motion for the averaged fields aL = (cil)^ 
<^R = (ci^r) read 



d 1 hv ■ 

— aL = - {-ex{t)aL + gaR) - -aL - \^b'^{t) 
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where 6^(t) = e~'^ ^ h^^ describes the electromagnetic 
drive along the feed line with amplitude 6*^ and frequency 
ujL' Here we used a rotating frame with laser detuning 
from resonance A/, = uol — ^q- The transmission T(t) = 




-2 2 

(a) Mean position x^/g (b) 



-4-2024 

Mean position x^/g 



FIG. 3. Transmission spectrum for the setup depicted in 
Fig.^for resonators' decay rate i^ — O.lg: density plot for the 
time-averaged transmission depending on mean mechanical 
displacement xo = exo, and frequency detuning Al = ujl—ooq 
of the feed line's microwave drive at ujl {ojq denotes the left 
mode's bare frequency for x = 0). (a) Without mechanical 
drive (x(t) = xo); the spectrum is given by the resonance 
frequency depicted in Fig. fill, (b) For mechanical driving 
(x(t) = Acos(Qt) + xo) with amplitude A — eA — 1.5Q and 
frequency Q = 3^; mechanical sidebands displaced by Q ap- 
pear and intersect the original photon branches where, due 
to mechanically driven Rabi dynamics, high transmission and 
additional anticrossings arise. The gap 2gi is determined by 
the Bessel function Ji according to 2gi — 2gJi[A/Q). 



K,{aj^{t)aR{t))/ (6*"') can be expressed as 



T{t) = n^ 



[ G{t,t')e 

-^ — CO 



-jAit'-(«/2)(t-i') ^^/ 



(6) 



where the phase comprises the feed line's drive and res- 
onators' decay, while the Green's function G(t, t') de- 
scribes the amplitude for a photon to enter the left res- 
onator's mode aL at time t' and to be found in the right 
one aR later at time t. 

We take into account two scenarios: first, the beam is 
at rest given a constant displacement x{t) = xq; second, 
the beam is mechanically driven to oscillate with ampli- 
tude A and frequency Vt around the mean position xq, 
x{t) = Xo + Acos(l]t). Fig.|3^ shows numerical results of 
the transmission spectrum without mechanical driving. 
The spectrum corresponds to the system's resonance fre- 
quency depicted in Fig. [T]i where the resonance width is 
set by the resonators' decay rate n. In contrast. Fig. [SfD 
shows the transmission including mechanical driving with 
Vt = 3^, i.e. g c:^ Vt being characteristic for coupled mi- 
crowave optomechanics. 

To understand the main features of Fig. ^p we note 
that, in general, two processes are involved to observe 
transmission, see ([6|: first, the left resonator aL must be 
excited by the electromagnetic drive 1)^^(1) = e~*^^^6*^; 
second, the internal dynamics must be able to transfer 
photons from aL to aR. From ([5| the solution G{t^t') 
can be found to be 



where (j){t') = {A/Q)^vii{Vtt') and aR{t^t') is a solution 
to the driven two state problem 
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with t >t^ and initial condition aL{t' ^ t') = 1, aR(t' ^t') = 
0. Note that we expressed displacement in terms of 
frequency; A = eA^ xq = cxq. For A 7^ 0, in ad- 
dition to the electromagnetic drive (see e~*^^ in (te|), 
the mechanical driving can excite aL in terms of mul- 
tiples of the mechanical frequency mQ. This mechani- 
cal excitation is described by the phase factor e~'^^^* ^ = 
Xlm Jm{A/^)e~^^^^' in (M) and leads to mechanical side- 
bands in the spectrum [cf. Fig. [sb]. Note that the in- 
dividual process mQ is described by a Bessel function 
Jm{A/rt) and can be tuned by the driving strength. Be- 
yond the modified excitation, the driving significantly 
changes the internal dynamics of the microwave fields, 
see Eq. Q. In particular the mechanical motion can 
initiate mechanically driven Rabi dynamics exchanging 
photons between aL and aR that leads to high transmis- 
sion if the mechanical drive at Q is in resonance with the 
modes' frequency difference. For sufficiently strong driv- 
ing, the mechanically assisted process leads to additional 
anticrossings in the spectrum resembling Autler-Townes 
splittings known from quantum optics (see marker in 
Fig. ^) . From Eq. (|8| we find that the spacing of this 
first additional splitting scales according to 2gJi{A/rt) 
and can likewise be tuned by the mechanical driving 
strength. All this illustrates how, due to ^ c^ Q^ the 
microwave field can extensively be manipulated by me- 
chanical motion in terms of mechanically driven coherent 
photon dynamics. 

Coupling to the square of displacement. - We 
present a modified scheme comprising coupled microwave 
resonators that allows to couple the photon number to 
the square of mechanical displacement [Fig. [4k-b]. In 
contrast to the setup in Fig. [T] here the nanomechanical 
beam is placed in the region between the two resonators, 
such that its motion affects both simultaneously [Fig.[4]3]. 
While for a given displacement x, the line capacitance c 
of the first resonator is increased, the one of the second 
wave guide is decreased and vice versa. According to 
our previous results, using the notation from above, the 
Hamiltonian for this setup reads 



H = h {ujQ — ex) a\aL + h {ujq 
+ hg ia^aR + a^^aLJ , 



I a^aR 



(9) 



G{t,t') = aR{t,t') 



-i4>(t') 



(7) 



where g — ujQCgdg/cd^ see (|3| and Q. 

Fig. I4J3 illustrates the system's resonance frequency 
cj^(x) = cjo ± \/g'^ + (ex)^ as function of displacement. 
Naturally all the characteristics of coupled multimode 
optomechanics in the microwave regime, that have been 
discussed above, apply. In particular the hyperbola- 
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FIG. 4. Schematic device geometry for two microwave res- 
onators aL , an with a nanomechanical device coupled to both 
of them, (a) Two striphne resonators (each of length d) are 
coupled to external feed and transmission lines (green). The 
central conductors of az,, clr (red) are capacitively coupled 
in a small region of length dg where the wave guides adjoin. 
(b) Between the two resonators a small mechanical beam, 
connected to ground (blue), is placed. Its displacement x af- 
fects the line capacitance of both, aL and an. (c) System's 
resonance frequency as function of displacement: the beam's 
displacement linearly changes the bare modes' frequency of 
aL and an (dashed). Due to the coupling g between the res- 
onators, there is an avoided crossing 2g in the eigenfrequen- 
cies uj± (blue), (d) Analogous optical setup with a movable 
dielectric membrane placed in the middle of a cavity j2j . (e) 
Schematic realization with two microwave LC circuits where 
a central plate is grounded and resonates against two others 
that build the LC circuits. In the notation of Fig. [2] we get 
for the coupling frequency g = uolcG/C with ulc — l/\/LC . 



shaped avoided level crossing allows to realize Landau- 
Zener transitions and the dynamics of Landau- Zener- 
Stueckelberg oscillations in the light field of a microwave 
setup. At the extrema of ijO±{x) the Hamiltonian allows 
an exclusive coupling to the square of mechanical dis- 
placement x^ . Thus, in the following, we investigate the 
prospects to perform QND Fock state detection using 
this microwave setting. 

Fock state detection. - In principle, an exclusive cou- 
pling of the photon number to x'^ allows to perform 
QND Fock state detection and to observe quantum jumps 
of a mechanical resonator [lOl [Tl]. Indeed the Hamil- 
tonian ([9]) corresponds to the one found for an opti- 
cal setup (see Fig. [iji), that generated a lot of inter- 
est in this regard [2, 12 . We focus on one microwave 
mode with annihilation operator a and expand the res- 
onance frequency uj-^{x) around xq. For xq = the 
linear contribution vanishes and the Hamiltonian reads 



H = h (a;+(0) + lujl{0)x%[b'^ + b]A a'^a + hnb'^b, where 

we considered the phonon number operator n = b^b and 
quantized (x — xq) using the mechanical beam's displace- 
ment operator x = Xzp (^^ + ^) . The zero-point displace- 
ment Xzp = ^/h/2mQ is determined by the mechanical 
mass m and frequency Q. Applying rotating wave ap- 
proximation (RWA), (b^ -\- b) c^ 2n-\-l^ we immediately 
see that [H^ n] = 0. For a potential experiment we con- 
sider a scheme in analogy to the one proposed for the 
optical setup [2 . The mechanics is cooled to the quan- 
tum mechanical ground state [23l [25] . After switching 
off the cooling, the phonon number n is measured via 
the frequency of the microwave mode. To detect a quan- 
tum jump from n = to n = 1, the frequency shift per 
phonon Au = uj_^{0)x1p, where uj_^{0) = e^/g, must be 
resolved within the lifetime of the phonon ground state 
r^^^. Given the imprecision of the frequency measure- 
ment in terms of the angular frequency noise power spec- 
tral density S^joj (in units s~^/Hz), the signal-to-noise 
ratio reads S = {Auj^t^^^ / S^joj [2^. In contrast to the 
optical regime, where shot-noise-limited frequency mea- 
surements are routinely achieved, microwave setups in 
general suffer from amplifier noise adding Uadd quanta 
of noise beyond the shot-noise limit in a Pound-Drever- 
Hall scheme, such that Soju = ij^add + l/2)H?fujOc/^^Pin 
(see [30]). Pin denotes the incident power and ujc is the 
resonance frequency of the cavity. While commercially 
available systems add a significant amount of noise, a 
new Josephson parametric amplifier achieved Uadd < 1/2 
^31] . This technique has already been used for displace- 
ment measurements in a microwave optomechanical sys- 
tem with Uadd = 1.3 [26]. Essentially, the total lifetime 
^(0) _ l/(r^"^+r^^^+r^Tj) is set by the thermal lifetime 
tt = hQ/ksT via the mechanical quality factor Q and 
the chip temperature T. Additional contributions due to 
the RWA (trwa) and imprecise positioning xq 7^ (run) 
will be determined via Fermi's golden rule rates (see [2]). 

For microwave setups using a small nanomechanical 
beam manufactured close to the central conductor of a 
stripline resonator, achieving optomechanical couplings 
of e c:^ 1 - lOOkHz/nm [HI [13 [25] , the frequency shift 
per phonon Auj ex e^ turns out to be extremely small 
making Fock state detection impossible. A new on-chip 
microwave system however, consisting of an LC circuit 
where the plates of a parallel-plate condensator mechan- 
ically resonate, achieves e = 65MHz/nm [32 . Our pro- 
posal transfers to this scheme by stacking three such 
plates, see Fig. |4^. For experimentally realistic parame- 
ters [32] , a calculation of S yields that a setup with this 
optomechanical coupling would allow to detect an indi- 
vidual quantum jump from the mechanical ground state 
to the first exited state, see Tab.[l| Note that, in contrast 
to the setup discussed in [2 , the parameters here are al- 
ready in the small-splitting regime g < Q^ and the details 
of Fock state detection in that regime may require further 



e [MHz/nm] 


m [pg] 


Q/27r [MHz] 


Q/10' 


K/n 


g/27r [MHz] 


P^n [PW] 


xo [pm] 
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T [mK] 
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65 


10 


11 


3.5 


1/70 
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20 
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70 


10 
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1/100 


0.5 


50 


0.5 
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20 
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TABLE I. Two sets of experimental parameters that would allow to observe an individual quantum jump from the mechanical 
ground state to the first excited state, with a signal-to-noise ratio S > 1. Further parameter a;c/27r = 5 GHz. 



analysis. Finally, we point out that such a setup, even for 
S < 1, would allow to measure "phonon shot noise", i.e. 
quantum energy fluctuations around an average phonon 
number, of a mechanically driven, ground-state-cooled 
mechanical oscillator [33 . 

Conclusion. - To conclude, we introduced and an- 
alyzed theoretically coupled multimode optomechanical 
systems for the microwave regime. In contrast to the op- 
tical domain, these systems possess coupling frequencies 
between the electromagnetic modes that are naturally in 
the range of typical mechanical frequencies {g c^ Q). By 
calculating the transmission spectrum, we demonstrated 
how this allows to manipulate the microwave field dy- 
namics in terms of mechanical driving. In principle g c^ Q 
enables to realize all kinds of driven two- and multi-level 
dynamics known from quantum optics in the microwave 
light field. Our discussion mostly focussed on classical 
mechanical motion. However, for mechanical oscillators 
in the quantum regime, coupled multimode systems with 
g c^ Q. will be particularly interesting. For instance it 
might be possible to realize hybridized states which are 
superpositions of states with a photon being in differ- 
ent modes and phonons in the mechanics. We further- 
more proposed a multimode setup that allows to couple 
the microwave photon number to the square of mechani- 
cal displacement and enables QND Fock state detection. 
For experimentally realistic parameters we predicted the 
possibility to detect an individual quantum jump from 
the mechanical ground state to the first excited state. 
The same scheme also allows to measure phonon shot 
noise. Both experiments would constitute a major break- 
through. 
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